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ABSTRACT 

We make use of our 'minimal' cold interstellar medium (ISM) emission line model that predicts the molec- 
ular and atomic line emission per unit dense, star-forming gas mass (Geach & Papadopoulos 2012; Paper 
I) to examine the utility of key line ratios in surveys of the so-called star formation 'mode' as traced by 
£sF=-^dense(H2)/M to tai(H2). We argue that £gp and its proxies provide very sensitive, extinction-free dis- 
criminators of rapid starburst/merger-driven versus secular quiescent/disk-like stellar mass assembly, with the 
most promising diagnostic to be applied in the near-future being CO 7(4— >3) I [C l]( 3 Pi-^ 3 Po) ■ These lines 
are accessible across nearly the full range < z < 2 (thus covering the bulk of galaxy evolution) with the 
Atacama Large Millimeter Array. In addition to their diagnostic power, another advantage of this combination 
is the similar observed frequencies (A^o ~ 30 GHz) of the lines, resulting in nearly spatially-matched beams 
for a fixed aperture, thus mitigating the effects of resolution/morphology bias in the interpretation of galaxy- 
averaged line ratios. Finally we discuss the capability of deep blind redshift surveys with the high frequency 
component of the Square Kilometer Arrray (SKA) in discovering H2-rich galaxies with very low £sf values. 
These could be the progenitors of starburst galaxies seen prior to the onset of star formation; such galaxies 
could be a class of extreme outliers from local (gas surface density )-(star formation rate) scaling laws, which 
would exclude them from current star formation or stellar mass selected samples. Our conservative model 
suggests that SKA could detect such systems residing at z ~ 3 at a rate of 20-200 per hour. 

Subject headings: galaxies: ISM — galaxies: starburst - 
ISM: molecules: CO, HCN 



galaxies: evolution — cosmology: observations ■ 



1. INTRODUCTION 

The prevailing 'mode' of star formation in the distant Uni- 
verse, namely merger-driven starbursts versus isolated star- 
forming gas-rich disks - with their H2 gas as traced by its 
carbon monoxide (CO) line emission - is now a topic of in- 
tense interest after the discovery of numerous galaxies at high 
redshifts with gas fractions above 50% that seem to follow the 
latter mode (e.g. Daddi et al. 2010; Tacconi et al. 2010). Nev- 
ertheless, the indicators used in such studies, namely the near 
Galactic value of the X co factor, and low-excitation low-/ CO 
spectral line energy distributions (SLEDs) are strongly degen- 
erate when it comes to the star formation mode (Papadopoulos 
et al. 2012) and thus cannot be used to securely argue for a bi- 
modal distribution of star formation modes. Moreover using 
ratios of single CO line luminosities (used as measures of to- 
tal molecular gas mass) to infrared luminosity (used as star 
formation rate proxy) as mode indicators is fraught with un- 
certainty because the value of the X co remains both highly un- 
certain and dependant on galaxy type, molecular gas surface 
density and kinematic state, as well as average star formation 
activity (see Narayanan et al. 201 1 for a recent exposition). 

In the local Universe the star formation 'mode' seems 
to be very well (i.e. uniquely) indicated by the fraction of 
the total gas reservoir residing in the dense phase: £gp = 
-^densc(H2)/A/totai(H2), a fact that can be easily exploited 
in high redshift galaxy surveys. In this article, we explore 
practical means of determining the star formation mode of 
galaxies, using low-/ CO or [C l]( 3 Pi— > 3 Po) and hydrogen 
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cyanide (HCN) line surveys as the most effective tool for de- 
termining £sF=M denS e(H 2 )/M to t a i(H2) as a function of red- 
shift. Selecting galaxies through £gp is independent of the 
various selection criteria used in traditional galaxy surveys, 
and is unaffected by degeneracies of CO SLEDs in determin- 
ing £sf (Papadopoulos et al. 2012). 

For the purposes of this work we seek good proxies of the 
dense gas mass fraction and investigate their observational 
prospects in the distant Universe, especially in survey mode, 
but will not delve into the details of how such proxies can be 
translated into actual £sf values. Multi-/ and multi-species 
molecular line observations are needed for this (e.g. Greve et 
al. 2009), which will soon become routine for the Atacama 
Large Millimeter Array (ALMA). We note that the prize from 
actual estimates of £§f values and direct measurements of tur- 
bulent line widths in the molecular gas of distant star-forming 
galaxies with sensitive imaging interferometers like ALMA 
(see Swinbank et al. 201 1 for an example for a gas-rich star- 
burst at z ~ 2.3) will be to properly investigate proposed uni- 
versal turbulence -regulated star formation theory for galaxies 
(Krumholz & McKee 2005). 

Finally we outline some of the unique discovery possibil- 
ities that very deep low-/ CO ( J up < 3) surveys with the 
Square Kilometer Array (SKA) can offer when it comes to 
gas-rich systems with very low £sf values that can only be 
discovered with blind surveys. Such galaxies will have negli- 
gible HCN/(1— >0) and CO J up > 4 line emission, and would 
be extreme outliers of the Schmidt-Kennicutt relation (Kenni- 
cutt 1998), unselected by current high-redshift galaxy surveys 
as these are based, directly or indirectly, on the star forma- 
tion rate (SFR, e.g. optical/near-IR lines from H II regions, 
(sub)mm dust thermal continuum, and cm non-thermal con- 
tinuum). This population would be valuable for extending 
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our understanding of the gas properties and early star forma- 
tion histories of massive galaxies. Where relevant, we assume 
a flat cosmological model, with O m = 0.3, fl\ = 0.7 and 
H Q = 70kms- 1 Mp C - 1 . 

2. PRACTICAL INVESTIGATIONS OF THE STAR FORMATION 
MODE OF GALAXIES 

2.1. Theoretical basis 

Paper I presented a model for the abundance of molecular 
and atomic line-emitting galaxies seen across cosmic time, 
using a 'minimal' model for the cold molecular interstellar 
medium (ISM). We envision that future surveys will assem- 
ble large samples of high-z systems with bright molecular 
line mission, be it in a totally blind sense, or by (sub)mm/cm 
spectroscopic follow-up of a SFR (most likely proxied by the 
submm continuum) or stellar mass limited sample (achievable 
with deep optical/near-infrared surveys). A survey of the star 
formation mode would be a powerful probe of galaxy evo- 
lution, since it would perform the dual task of measuring the 
evolution of the total gas fraction of star-forming galaxies, the 
nature of the star formation (burst or secular), and a more ac- 
curate prediction of the rate of consumption of the cold gas 
reservoirs. 

A prominent bi-modality of the £sf values in merger-driven 
versus isolated disk star formation is actually expected as 
a general characteristic of supersonic turbulence because of 
the much larger (^5-30 x) turbulent line widths observed in 
CO line emission in mergers/ultraluminous infrared galaxies 
(ULIRGs) than in isolated spiral disks (Downes & Solomon 
1998). This difference is expected to be a general charac- 
teristic of merger versus isolated disk systems (Narayanan et 
al. 2011). The resulting difference in the corresponding 1- 
dimensional average Mach numbers Aiulirgs ~ (5 — 30) x 
■M disks (where M = a v /c s ) leads to much higher dense 
gas mass fractions in the ISM of mergers than in isolated 
disks given that the probability distribution function of the 
molecular gas density in supersonic turbulence is well ap- 
proximated by a log-normal distribution with a dispersion of 
<r p = y/ (In (1 + 3M 2 /4)) (Padoan & Nordlund 2002). 

The very enhanced dense gas mass fractions molecular gas 
mass fraction expected to be in structures with overdensities 
x > x (where x = n/(n), and n is the gas density) is then 
given by 
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The typically large velocity dispersions of (7^(30- 
140) km s -1 in the molecular gas disks of the merger systems 
found in ULIRGs (and recently even in a molecular gas disk 
fueling star formation in an spectacular starburst at z = 2.3, 
see Swinbank et al. 2011), with corresponding large Ai will 
dramatically extend the gas density probability distribution 
function towards high values where much of the molecular 
gas mass will lie, quite unlike the less turbulent molecular 
gas in isolated disks (see Figure 1). Thus for giant molecular 
clouds (GMCs) in typical spiral disks only ^3% of their mass 
is expected to be at overdensities x > 500 (5x 10 4 cm~ 3 for 
a typical GMC with (n) = 100 cm~ 3 ) while for the highly 
turbulent gas in mergers this is ~45%, and can be even higher 
(e.g. Greve et al. 2009). 
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Figure 1. The gas mass fraction expected at overdensities of x > xq (where 
x = n/(n), see Equation 1) in Giant Molecular Clouds (GMCs) for a 
given 1-dimensional average Mach number M = cr v /c s . The large boost 
of M(x > £o)/A^total that occurs in GMCs resident in merger/ULIRGs 
(M ~ (30 — 300)) with respect to those in spiral disks (A4 ~ 10) is evi- 
dent. 



2.2. Practical tracers of the star formation mode 

In the following we consider various combinations of 
molecular and atomic lines as proxies of £sf- A summary 
from available line data for mostly local galaxies can be found 
in Figure 2 where we plot normalized line ratios for a sample 
of observations of (i) quiescent star forming disks and dark 
clouds and (ii) galactic nuclei, mergers and starburst galax- 
ies, compiled from the literature (Gao & Solomon 2004, Is- 
rael et al. 1998, 2001, 2003, Barvainis et al. 1997, White et 
al. 1994, Petitpas & Wilson 1998). We define the normalized 
ratios as r a /b/( r a/b)^ where a and b correspond to (for exam- 
ple) the line brightness luminosities of the dense and total gas 
tracers HCN 7(1^0) and CO/(1->0) with (r a / b ) represent- 
ing the average over the entire sample. 

It is clear from Figure 2 that carefully-picked line diag- 
nostics can sensitively distinguish between 'disk' and 'burst' 
mode, with r a / b /(r a / b ) « 1 dividing the two classes. In- 
dividual £ S f tracers (e.g. HCN7(1— >0) and CO 7(1^0) line 
ratios) can be even more sensitive than what the distribu- 
tion in Figure 2 indicates, with factors of ^10 difference be- 
tween 'disk' and 'burst' mode (see §2.2.1). Moreover, the 
sample presented in Figure 2 is dominated by local systems 
while high redshift observations present more of a challenge 
as global line measurements may represent wide spatial aver- 
ages (smearing out the line ratio for composite systems), but 
also because the observed frequencies of lines available in the 
local Universe will be redshifted out of the current range of 
sensitive facilities for certain important cosmological spans 
of time. 

We now will now discuss several of these tracers in more 
detail, highlighting their relative merits with a view to iden- 
tifying diagnostic tools that will be suitable for practical sur- 
veys of the star formation mode of galaxies over large spans 
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of cosmic time. 
2.2.1. HCN7(1^0)/COy(1^0), HCN/(2->l)/CO/(2->-l) 

This is the cleanest £sf indicator available, with local stud- 
ies showing a clear bi-modality for the ISM in merger-driven 
starbursts versus isolated disks with lower SFRs in the lo- 
cal Universe (Gao & Solomon 2004), with HCN tracing the 
dense, star-forming gas mass, and low-/ CO tracking the to- 
tal gas reservoir. Arguably this line combination is the most 
direct tracer of £sf; extensive observations of local galaxies 
show this line ratio to be an excellent way to discriminate be- 
tween different modes, with high values of ^hcn/co ~ 0.2- 
0.3 found exclusively in compact, merger-driven, extreme 
starbursts while ?*hcn/CO ~ 0.01-0.03 for isolated spiral 
disks with extended star formation (e.g. Solomon et al. 1992; 
Gao & Solomon 2004). 

Observing HCN 7(l->0) O = 88.632 GHz) and 
CO/(1^0)O = 115.271GHz) has the disadvantage 
that, although this line combination will be accessible with 
ALMA for local systems (important for an even more 
robust calibration in terms of £sf)> by z ~ 1 CO 7(1— )-0)is 
redshifted out of the ALMA bands and remains inaccessible 
until z ~ 1.5 in the cm bands (and accessible with other 
cm-wave facilities). In addition, due to the broad telluric 
absorption feature at v ~ 60 GHz, there will be a broad 
HCN 7(1— ^0) and CO 7(1-^0) 'desert' in the z < 1 interval, 
limiting the use of this tracer for intermediate redshift work. 

Due to the aforementioned shortcomings, this particular in- 
dicator will only be effective for probing very local galaxies 
rather than those at cosmologically significant redshifts of z > 
2, and thus cannot satisfactorily measure £sf over the peak of 
the global star formation rate density in a well-sampled, to- 
mographic manner. However, given the small range of both 
CO and HCN r 2 \ ~ 0.6-1, the J = 2 -> 1 transition of these 
two molecules can also be used to the same effect. This allows 
ALMA to survey r HCN / co up to z < 1 as the rest frequen- 
cies of CO 7(2->l) and HCN 7(2->l) are at 230.538 GHz and 
177.264GHz respectively. HCN7(2— >l)can be accessed by 
ALMA out to z ~ 1, and HCN 7(3^2) out to z ~ 2. 
Given the order of magnitude differences expected for the 
r HCN/co"tyP e of ratios between mergers and isolated disks 
it is very likely that even HCN 7(3—5-2) and ratios such as 
HCN 7(3^2) / CO 7(l->0) or HCN 7(3^2) / CO 7(2->l) will 
remain practical proxies of £gp out to z < 2 using ALMA. 
Moreover, given the tremendous sensitivity boost at (sub)mm 
wavelengths that ALMA represents, the necessary calibration 
of such 'hybrid' HCN/CO ratios in terms of £sf will be trivial 
for large numbers of star-forming galaxies in the local Uni- 
verse. 

2.2.2. HCN7(l->0)/[Cl]( 3 Pi-> 3 P ), 
HCN7(2-»1)/[C l]( 3 Pi^ 3 P ) 

This is a potentially powerful tracer of £sf, with - as men- 
tioned above - the HCN7(1^0)and HCN 7(2-5-1) line lu- 
minosities proportional to the dense star-forming gas mass, 
and [C l]( 3 Pi— > 3 Po) scaling with total molecular gas mass, 
but with many advantages over low-7 CO lines (see Pa- 
per I). The higher rest frequency of [C l]( 3 Pi— > 3 Po) (^o = 
492 GHz), compared to HCN 7(l->0) means that the corre- 
sponding [Cl]( P\— > 3 Pq) observations would be conducted 
with ALMA Bands 6 and lower. 

Local data support the notion of the ratio of 
HCN 7(1^0) and [C l]( 3 Pi^ 3 P ) lines as a £ SF 



proxy, with rHCN/[ci] ^0.55-2.2 observed in ULIRGs 
( [C l]( 3 Pi— > 3 Po) data presented in Papadopoulos & Greve 
2004; Greve et al. 2009, and HCN 7(1^0) data from Gao & 
Solomon 2004). In comparison, for disk-dominated quiescent 
galaxies, /"HCN/tCil ~ 0.06 — 0.10, although it can rise 
up to ^0.5 for vigorously star-forming disks like that of 
NGC1068 (using HCN 7(1^0) data from Gao & Solomon 
[2004] and [C l]( 3 Pi^ 3 P ) data from Israel et al. [2009]). 
Further benchmarking work on this ratio as a practical 
mode indicator must be done locally (now possible from 
high altitude dry sites like the Chajnantor plateau) before 
applying it in the distant Universe. A larger body of multi-7 
HCN line data along with [C l]( 3 Pi-> 3 Po) measurements 
and models of the [C/H2] abundanc^ are necessary to 
both benchmark r HCN /[ci]-typ e ratios as £sf indicators as 
well as help disentangle the effects of high SFRs (and thus 
[C/H2]-boosting cosmic rays) on this valuable mode indicator 
(see also §2.2.4). 

2.2.3. 12 CO/ 13 CO ratios for J up < 2 

The Pi2/i 3 (l -> 0, 2 ->• 1) = 13 CO/ 12 CO ratios are not 
proxies of £sf but rather of the turbulent velocity fields of 
giant molecular clouds in star-forming galaxies. Numerous 
studies have shown that ^12/13(1 — > 0, 2 — > 1) ~ 5 — 15 
in disks, with low values found in Milky Way type systems 
and larger ones in vigorously star-forming disks (e.g. Casoli 
et al. 1992; Aalto et al. 1995; Papadopoulos & Seaquist 1999). 
Ratios ^12/13(1 — > 0, 2 — !• 1) > 20 are found exclusively in 
merger systems and are thought to be the result of the highly 
turbulent ISM and the enhanced gas temperatures produced 
by the dissipation of their supersonic turbulence (Aalto et al. 
1995; Greve et al. 2009 and references therein). 

The advantage of using isotopologue line ratios as prox- 
ies of merger versus isolated-disk star formation modes, even 
if such ratios are not proxies of £gp, is obvious: the con- 
venience yielded by their frequency proximity. Wide band- 
width receivers allow the possibility of deep surveys with si- 
multaneous CO and 13 CO line observations. Objects with 
^12/13(1 — > 0, 2 — s- 1) > 20 would be considered as har- 
bouring merger-driven star formation, while all the rest would 
be in the disk mode. This method can safely be used up to 
J up = 2, partly because of lack of local 13 CO 7(3— > 2) or 
higher-7 line data (the 13 CO 7(3— >2) line lies precariously 
close to an atmospheric absorption feature) but also because 
higher-7 CO lines will progressively cease to be good tracers 
of the average velocity fields of molecular clouds and begin 
tracing their self-gravitating, dense regions (Ossenkopf 2002). 
Nevertheless recent measurements of a large Pi2/i3(7 = 
3 -> 2) in a distant ULIRG-hosting QSO (Henkel et al. 2010) 
indicate that such a diagnostic may remain a good star forma- 
tion mode indicator up to J up = 3. 

Finally we note that a competing explanation for the very 
high i?i2/i3(l 0, 2 1) ratios observed in all ULIRGs 
exists, attributing them to much higher [ 12 C/ 13 C] abundances 
in such systems (Casoli et al. 1992; Henkel & Mauersberger 
1993 and references therein), possibly caused by a different 
stellar initial mass function in merger/starburst systems. If 
this is true it would make P12/13 a less 'clean' star formation 
mode indicator unless a strong link can be proven between 

4 It must be stressed that abundance uncertainties affect all H2 mass mea- 
surements utilizing optically thin lines (e.g. 13 CO). Among those atomic car- 
bon is the species with the simplest, mostly cosmic ray controlled chemistry. 
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Figure 2. Comparison of the £gp proxies 7"hcn/CO> r HCN/[Cl] an d 
^00(4— 3)/ [CI] compiled for galaxy disks and clouds (representing our 'disk 
mode') and mergers and nuclear regions (representing our 'burst mode'); data 
taken from Gao & Solomon (2004), Israel et al. (1998, 2001, 2003), Barvai- 
nis et al. (1997), White et al. (1994), Petitpas & Wilson (1998). To present 
the tracers on the same scale, we have normalized each proxy to the average 
value, r" a /b/{ r a/b}> where a and b correspond to (e.g.) HCN7(1— >0) and 
CO 7(1— >0) . Although there is some overlap, due to the J"co(4— 3)/[CI] 
r HCN / [CI] ratios (see discussion in §2.2.2 and §2.2.4) there is a clear distinc- 
tion between the two classes of systems irrespective of the particular indicator 
used, with r a/b / (r a/b ) < 1 for quiescent disks and r a/b /{r a/b ) > 1 for 
merger/nuclear systems. 

such isotopic enhancements and merger-driven star-forming 
systems. 

2.2.4. COy(4^3)/[Cl]( 3 Pi^ 3 P ) 

The rco(4-3)/[Ci] rat i° has the key advantage that the lines 
used have similar rest frequencies, (Afrj « 30 GHz, and thus 
similar beams for a given aperture), while both remain acces- 
sible to ALMA across a large span of cosmic time. It is as- 
sumed that CO 7(4— >3) traces mostly dense star-forming gas, 
an expectation borne out both by observations (e.g. Petitpas & 
Wilson 1988; Nieten et al. 1999; Papadopoulos et al. 2012) as 
well as models of emergent CO line emission from turbulent 
molecular clouds (Ossenkopf 2002). Finally, both lines ben- 
efit from a positive ^-correction with respect to lower-/ lines 
making them suitable for (H2 mass)-sensitive high-z work. 

Unfortunately, their relatively high frequencies and the 
lack of multi-beam receivers at such frequencies has pre- 
vented the assembly of a large body of CO 7(4— !>3) and 
[C l]( 3 Pi— > 3 Po)line data in the extragalactic domain. Nev- 
ertheless, enough data is assembled to examine if a signifi- 
cant bi-modality of ?"co(4-3)/[ci] indeed exists for starburst 
environments dominated by large dense gas reservoirs (typi- 
cal of merger-driven systems) versus low-SFR environments 
in disks. In this comparison we also consider CO 7(4— > 3) and 
[Cl]( 3 Pi-> 3 P )data obtained for galactic nuclei (e.g. NGC 
253) including that of the Galactic Center as representative 
of the average molecular ISM properties of the merger-driven 
starbursts in ULIRGs (the high turbulence yielding also very 



high Afdense/^ftotai values in galactic centers). For such en- 
vironments we obtain (?'co(4-3)/[ci]) — 4.55 ±1.5 (data 
from Israel et al. 1995, 1998, 2001, 2003; Petitpas & Wil- 
son 1998; Barvainis et al. 1997; Papadopoulos & Greve 2004; 
Papadopoulos et al. 2004) where the dispersion represents 
that of the measured values. These include the Galactic Cen- 
ter, the nuclei of various nearby galaxies (NGC 253, Maf- 
fei2, M83), a distant strongly-lensed QSO (The Cloverleaf 
at z = 2.56) and a local ULIRG/QSO (Mrk231). For disk- 
dominated star-forming environments, the available literature 
data yield (f*co(4-3)/[ci]) ~ 0.45-1.3 with the lowest value 
(—0.43) found for the outer Galaxy and a quiescent cloud 
in M31 (<0.45). Somewhat higher values are obtained for 
the inner Galaxy (—0.83) and star-forming spirals such as 
NGC 4826 and NGC 3079 (-1.3, and see Figure 2). 

Thus it is seems that the ?"co(4->3)/[ci] rat io can be a practi- 
cal star formation mode indicator in the ALMA regime, with 
differences up to factors of —10 between mergers and disk- 
dominated systems. We note however that the intense cos- 
mic ray fluxes expected in starburst environments can signif- 
icantly boost the [C1/H2] abundance ratio (e.g. Papadopoulos 
et al. 2004 and references therein), reducing the ?"co(4-s.3)/[ci] 
ratio in extreme starbursts. This may be happening in the 
inner 500 pc of M82 where a rather low ?"co(4^3)/[ci] = 
1.5 is measured (White et al. 1994). Chemical models of 
[C1/H2] as a function of average cosmic ray energy densi- 
ties (cxSFR density) are necessary to account for such ef- 
fects as well as reducing the [C1/H2] abundance uncertain- 
ties (also needed for M to t a i(H2) estimates when using the 
the [C l]( 3 Pi— ^ 3 Prj) line luminosity). We expect that even 
crude models of a SFR-dependant (i.e. cosmic ray-dependant) 
[C1/H2] average abundance will yield ?"co(4^3)/[ci] (cosmic 
ray=Galactic) ratios (i.e. referenced to a common Galactic 
cosmic ray energy density) that will be even more sensitive 
to £sf (and thus the star formation mode) by much reduc- 
ing the little overlap that exists for this currently cosmic ray- 
uncorrected star formation mode proxy (and r HCN /[ CI ]) in 
Figure 2. 

2.3. C07(4— s>3)and [C l]( 3 Pi—!- 3 P ) line flux predictions 
and the design of a practical star formation mode survey 

In light of the distinct advantages that this diagnostic of- 
fers for high-z studies, in Figure 3 we plot the fluxes of 
the [Ci]( 3 PW 3 P )and CO 7(4^3) lines for a L m = 5 x 
10 U L@ galaxy at < z < 2 undergoing a burst (£sf = 0.6) 
and disk-mode star formation (£sf = 0.04). To calculate the 
line fluxes, we use our minimalist ISM model presented in 
Paper I, assuming the 'supervirial' and 'virial' versions of the 
model for the burst and disk mode respectively (see §3 of Pa- 
per I for further details). The [C l]( 3 Pi — > 3 7o) line flux is es- 
timated from the total gas mass following equation 1 1 of Pa- 
padopoulos, Thi & Viti (2004), adopting an excitation value 
of Qio = 0.35 (where Qiq — ATio/iVrci] is the ratio of the 
column density of gas in level 3 Pi to the ground state, ap- 
propriate for the lower temperature and excitation conditions 
when averaged over a global gas reservoir. The line ratios pre- 
dicted from our model are 7"co(4-s.3) /ci = 0.3 for £sf = 0.04 
and ^co(4-i-3)/ci — 2.1 for £sf = 0.6; in agreement with the 
normalized values presented in Figure 2. 

We compare the predicted line fluxes to the 5 a sensitivity 
(for channels of 300 km s -1 ) of full-power ALMA (assumed 
to be 50 antennas) in integrations of 1, 10 and 24 hours. This 
highlights the suitability of this line combination as a tool for 
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Figure3. Predicted CO 7(4— >3) and [C 1]( 3 P 1 ^- 3 P ) line fluxes for a 
galaxy with Lir = 5 x lO 11 !/© in the burst mode (£sf = 0.6) and disk 
mode (£sf = 0.04) seen at < z < 2. Our line flux estimates are based on 
our 'minimal' ISM model presented in Paper I. The benefit of using this line 
combination as a proxy for £gp is in the similarity of line frequencies, and so 
spatially matched beams for a given telescope, and the near-continuous ac- 
cessibility of the lines over cosmological distances from the ground (the blue 
shaded region shows the atmospheric transmission - on an arbitrary vertical 
scale - for a column of 1 mm of precipitable water vapour at Chajnantor, in- 
dicating the available windows for observations). Here we show the 5<r sen- 
sitivity of full power ALMA (50 antennas) achieved in 1, 10 and 24 hours, 
demonstrating the feasibility of performing a star formation mode survey for 
a luminosity-selected sample of galaxies over the majority of cosmic history, 
covering the peak epoch of star formation to the present day. 

surveys of the star formation mode, with both CO 7(4— >3) and 
[C l]( 3 Pi— ?> 3 Po) li nes detectable in realistic times for rela- 
tively normal galaxies over the full span of cosmological time 
pertinent to galaxy evolution in both disk and burst modes. 
The similarity of the line frequencies are key, because they 
provide nearly spatially matched beams (for a given dish, and 
in the case of interferometers, array configurationsQ), and this 
can be critical in mitigating the effect of any structural bias in 
the interpretation of such line ratios due to the 'mixing' of line 
emission from systems with different star formation modes 
that may exist in relative close proximity in merger configura- 
tions (e.g. Ivison et al. 2010). Moreover, in the cases of well- 
resolved individual systems the ?"co(4^3)/ci rat i° can g lVe 
a quick approximate 'reading' of the Ed C nso/S(H2) gas sur- 
face density ratio as a function of position within the galaxy. 
This level of detail, while certainly of interest to other types 
of studies concerning, for example, the detailed gas dynamics 
and its inward transport in galactic disks, will still be valu- 
able for uncovering 'milder' merger systems where the orig- 
inal disk is not destroyed but its Md nsc(H2)/Mt o tai(H2) is 
nevertheless enhanced by minor merger events. 

3. SEARCHING FOR H 2 -RICH GALAXIES WITH ULTRA-LOW STAR 
FORMATION RATES IN THE DISTANT UNIVERSE 

5 In Paper I and here, we have assumed the most compact array configura- 
tion will be used in survey mode. 



One of the most exciting aspects of deep blind molecular 
line surveys is their potential to discover gas-rich systems 
with very low levels of star formation; i.e. very low £sf- 
This could be a new class of objects that may be currently 
omitted by even the deepest galaxy surveys as these depend 
always on the SFR either directly (e.g. optical/near-infrared 
lines from H II regions) or indirectly (non-thermal cm contin- 
uum, (sub)mm dust continuum). The SKA and its pathfinders 
can search for the H I phase of such systems while deep low- 
/ CO and [C l]( 3 Pi— s- 3 P ) line observations (the latter with 
ALMA) can reveal their molecular gas reservoir. Of course 
the existence of the latter and its detectability via CO and 
[C l]( 3 Pi— > 3 Pq) lines assumes some significant metal enrich- 
ment of its ISM (so H2 can form within reasonable time and 
CO, [Cl]( 3 Pi— > 3 P )can be there to trace it) and thus some 
initial star formation episode. This need not contradict very 
low SFRs in the subsequent evolution of such systems, as a 
combination of structural properties (e.g. large thin disks), 
star formation feedback, and very different average ISM prop- 
erties (very low density molecular clouds with little molecular 
cooling) may act to suppress SFRs over long periods of time. 

This has already been shown possible, albeit for short pe- 
riods of time, for gas-rich and metal-poor disks with other- 
wise typical structural characteristics (Papadopoulos & Pelu- 
pessy 2010). During such periods these systems can undergo 
'bursts' of very inefficient SF (as well as 'bursts' of very ef- 
ficient star formation) that would place these systems as ex- 
treme deviants of the local Schmidt-Kennicutt (S-K) SFR/gas 
density scalings. Regardless whether H^-rich galaxies with 
ultra-low SFRs are simply a) short (S-K)-deviant periods over 
much longer lifetimes of otherwise normal SF galaxies (e.g. 
'pre-starbursts': gas-rich systems seen just prior to a burst 
of intense star formation), or b) a completely new class of 
objects widely spread across cosmic epoch, deep low-/ CO 
and [C l]( 3 Pi— !> 3 Po) blind line surveys are the way to find 
them. The [C l]( 3 Pi— > 3 Po) line is particularly promising be- 
cause of its large positive ^-correction (with respect to low- 
/ CO lines) yet easy excitation, and its potential to trace the 
metal-poor ISM with much of CO dissociated by far-UV radi- 
ation (or simply never formed in sufficient abundances) leav- 
ing [Cl]( 3 Pi— > 3 P )and [Cll] as the dominant form of car- 
bon. On the other hand, the SKA will have the sensitiv- 
ity and field of view making it favourable for searching for 
low-/ CO at high-redshift compared to [C l]( 3 Pi-J> 3 P )in 
the ALMA bands. We note that the [C II] line can trace such 
ISM much more sensitively than [C l]( 3 Pi— >- 3 Po)but unfor- 
tunately it becomes accessible to sensitive ALMA bands only 
for z > 2 (Figure 2 in Paper I), and it cannot be easily trans- 
lated to total molecular gas mass like the latter. 

Here we consider a naive model, where it is assumed that 
every star-forming galaxy goes through a brief pre-starburst 
phase prior to the onset of star formation, such that the abun- 
dance of sources is simply proportional to number of galax- 
ies contributing to the infrared luminosity density used in our 
abundance model (Paper I): $psb = A$ir(^, L). We assume 
a pre-starburst 'duty cycle' of A = 0.1 — 1% and 'low' mode 
of £sf = 0.01, such that the majority of gas is in the cold, qui- 
escent phase. Clearly the SKA with its sensitivity, wide field 
of view and access to low-/ CO lines at z > 3 is the most 
suitable facility for a search of pre-starburst galaxies. We 
present the cumulative number counts of CO /(l— >0) emitters 
at fobs = 30 GHz (assuming an instantaneous bandwidth of 
4 GHz) and A = 0.1 — 1% in Figure 4. Operating at the 'opti- 
mal' flux limit that maximizes the detection rate for a fixed ob- 
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Figure 4. Predicted abundance of CO7(l—s>0) emitters at u a ^ s = 30 GHz 
representing 'pre-starburst' galaxies, i.e. gas rich galaxies with £sf = 0.01 
(i.e. very low values) at z 3. The abundance estimates are calculated using 
the the model presented in Paper I. Our model ties the number density of 
CO emitters to the evolution of the infrared luminosity density (Befhermin 
et al. 2011), relating the luminosity of the principle molecular and atomic 
lines using our minimal ISM model. To predict the abundance of these pre- 
starbursts, we assume a naive model where the space density of pre-starbursts 
is related to the infrared luminosity density like 'I'psb = A<E>ir(z,L). In 
this plot we show the CO 7(1— ^0) surface density for a range A = 0.1— 
1%. Assuming the latest design specifications for the SKA high frequency 
component, operating at the 'optimal' flux limit - where the counts become 
proportional to F~ 1 - we expect the blind discovery of 20-200 of these S-K 
outliers per hour. 

serving time (i.e. where the integral counts are proportional to 
F - 2 , see Blain et al. 2000; Carilli & Blain 2002), we predict 
that a blind survey with SKA could detect these S-K outliers 
at z s» 3 at rate of 20-200 galaxies per houiQ- This model is 
certainly too simplistic; cosmological simulations that include 
a prescription for the chemistry and the evolution of the cold 
molecular gas might be better placed to make predictions for 
the abundance and detectability of S-K outliers with ultra-low 
SFRs at high-z (Obreschkow et al. 2009; Obreschkow, Hey- 
wood & Rawlings, 201 1; Lagos et al. 201 1, 2012; J. E. Geach 
et al. 2012 in preparation). 

Finally we note that imaging of dense gas tracers such 
as HCN and carbon monosulfide (CS) line emission along 
with indicators of total molecular gas distribution such 
as CO/(1->0), CO/(2^1)and [C l]( 3 Pi^ 3 P ) lines is 
of paramount importance for extracting information on 
important evolutionary parameters of galaxies such as 
Miyn/-Mtotai(H2) (with more gas-rich systems expected at 
high redshifts irrespective of the star formation mode) and 
gas consumption timescales rgp = Md e nse(H2)/SFR. In- 
deed it is only for the dense gas phase that the latter have the 
intended meaning as the duration of the observed star forma- 
tion event, though in the literature the longer timescales of 
r ~ Af to tai(H2)/SFR are often used, a rather poor choice as 

6 assuming the latest design specifications for a high-frequency SKA: 
skatelescope.org, S. Rawlings, private communication, 2011 



not all of the molecular gas is relevant to star formation, only 
the dense phase. 

4. CONCLUSIONS AND SUMMARY 

In this work we use our minimal molecular and atomic line 
emission model presented in Paper I (where it was used to as- 
sess the practicality of blind galaxy surveys using ISM molec- 
ular and atomic lines), to discuss practical surveys of the star 
formation 'mode' of galaxies, i.e. isolated disks versus merg- 
ers, using the mass fraction £gp of the dense star-forming gas 
phase versus the total molecular gas mass residing in star- 
forming galaxies. We argue that theoretical expectations as 
well as local molecular line studies point to £sf as a very sen- 
sitive, extinction-free, measure of the star formation mode, 
independent of any galaxy morphological criteria. Moreover 
£sf does not suffer from well-known degeneracies of low-/ 
CO lines and the values of the so-called Xco factor which 
make them poor star formation mode indicators of galax- 
ies. Among the various proxies of £gp we discuss, includ- 
ing HCN/CO (J up < 2), HCN7(1^0)/[Cl]( 3 Fi^ 3 P ), 
ratios, we find that the CO7(4->3)/CO7(l->0) ratio cur- 
rently offers the most immediately practical path to star for- 
mation mode surveys using ALMA, with the similarity of 
the lines' frequencies offering observing efficiency, and the 
nearly matched beam size for a given telescope (the most suit- 
able being ALMA for galaxy evolution studies) being an im- 
portant advantage. Sensitive 12 CO/ 13 CO for J up < 2 line 
ratio surveys are another promising route maintaining the lat- 
ter advantage, though such ratios are not direct proxies of £gp 
rather than of the average turbulence levels in the ISM envi- 
ronments of disks (with low 12 CO/ 13 CO ratios) and mergers 
(with high 12 CO/ 13 CO ratios). 

We note that for uniquely parametrizing the £sf in gas- 
rich star-forming systems in the distant Universe, the specific 
high-density gas tracers used are not important. As sensitiv- 
ities dramatically improve in the era of ALMA, other tracers 
will become available; for example, CS lines will also be- 
come suitable, with the same merger-versus-disk bi-modality 
expected for the corresponding rcs/co an d res/ roil ratios. 
Of course before any foray into the high redshift Universe 
any proxy of £sf must be calibrated locally with large sam- 
ples of star-forming galaxies, a task easily performed with the 
vastly improved sensitivities of ALMA for all typical dense 
(and total) molecular gas mass tracers. 

Finally we comment on the capabilities of upcoming ob- 
serving facilities in searching for H2-rich galaxies with ultra- 
low star formation rates in the distant Universe. These would 
be extreme outliers of the local (SFR density)-(gas mass) 
Schmidt-Kennicutt (S-K) relations, and deep low-/ CO and 
[Cl]( 3 Pi-> 3 P )line observations with the SKA and ALMA 
are the only way to discover them. The [C II] line can also be 
used to trace such systems very effectively with ALMA, es- 
pecially if they are metal-poor, but only at very high redshifts 
where the line is redshifted into the sensitive ALMA bands 
(Paper I), while its luminosity cannot be straightforwardly as- 
signed to total molecular gas mass. Future cosmological sim- 
ulations of galaxy evolution that include their molecular gas 
phase can inform on such possibilities, but it must be empha- 
sized that such a putative population of (S-K) relation outlier 
gas-rich galaxies with ultra-low SFRs represents pure discov- 
ery space for deep molecular and atomic ISM line surveys us- 
ing ALMA, SKA and its precursors, and thus should certainly 
be undertaken. 
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